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ABSTRACT 

In this paper we describe the predictions for the smoothed weak lensing shear 
and aperture-mass of two simple analytical models of the density field: the minimal 
tree-model and the stellar model. Both models give identical results for the statistics 
of the 3-d density contrast smoothed over spherical cells and only differ by the de- 
tailed angular dependence of the many-body density correlations. We have shown in 
previous work that they also yield almost identical results for the pdf of the smoothed 
convergence, n s . We find that both models give rather close results for both the shear 
and the positive tail of the aperture-mass. However, we note that at small angular 
scales (9 S ^ 2') the tail of the pdf P(M ap ) for negative M ap shows a strong variation 
between the two models and the stellar model actually breaks down for 9 S 55 0.4' and 
M ap < 0. This shows that the statistics of the aperture- mass provides a very pre- 
cise probe of the detailed structure of the density field, as it is sensitive to both the 
amplitude and the detailed angular behaviour of the many-body correlations. On the 
other hand, the minimal tree-model shows good agreement with numerical simulations 
over all scales and redshifts of interest, while both models provide a good description 
of the pdf V(^is) of the smoothed shear components. Therefore, the shear and the 
aperture-mass provide robust and complimentary tools to measure the cosmological 
parameters as well as the detailed statistical properties of the density field. 
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1 INTRODUCTION 

Magnification and shearing in the images of high-redshift 
galaxies arise naturally due to gravitational lensing effects. 
These effects derive from the fluctuations in the gravita- 
tional potential related to the underlying density field. Sta- 
tistical analysis of such observed weak lensing data is there- 
fore very effective in probing the underlying density, which 
is assumed to be dominated by the dark matter content. 
Observational surveys (see, e.g., Bacon, Refregier & Ellis, 
2000, Hoekstra et al., 2002, Van Waerbeke et al., 2000, and 
Van Waerbeke et al., 2002) have consequently been partic- 
ularly fruitful in this regard and have enabled estimates for 
the cosmological parameters to be made. 

To model weak gravitational lensing, cosmological N- 
body simulations have been widely used in which the par- 
ticle and mass distributions have been set to reflect those 



expected in the real universe. The first numerical studies 
for gravitational lensing employing iV-body simulations used 
ray-tracing methods (see, e.g., Schneider & Weiss, 1988, 
Jarosszn'ski et al., 1990, Wambsganns, Cen & Ostriker, 
1998, Van Waerbeke, Bernardeau & Mellier, 1999, and Jain, 
Seljak & White, 2000) to follow the deflections of light from 
sources at high redshift in the simulations. More recently, 
Couchman, Barber & Thomas (1999) developed an algo- 
rithm to compute the full 3-dimensional shear matrices at 
locations along lines of sight throughout simulation volumes. 
Barber (2002) extended the method to combine the shear 
matrices in the appropriate fashion along the lines of sight 
to produce final Jacobian matrices from which the weak lens- 
ing statistics could be directly derived. 

On large angular scales, analytical computations for 
weak lensing statistics can readily be made, as this is the 
regime where perturbative calculations apply (e.g., Villum- 
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sen, 1996, Stebbins, 1996, Bernardeau et al., 1997, Jain & 
Seljak, 1997, Kaiser, 1998, Van Waerbeke, Bernardeau & 
Mellier, 1999, and Schneider et al., 1998). However, on small 
angular scales, especially relevant to observational surveys 
with small sky coverage, perturbative calculations are no 
longer valid and models to represent the gravitational clus- 
tering in the non-linear regime have been devised. 

To describe the non-linear evolution of the matter power 
spectrum, Hamilton et al. (1991) proposed a technique based 
on a non-local transformation. Their method was extended 
by Peacock & Dodds (1996) and included the conservation 
of mass and the rescaling of physical lengths in the different 
regimes, following the "stable-clustering" Ansatz of Peebles 
(1980). More recently, Peacock & Smith (2000) and Seljak 
(2000) developed the "Halo Model" for the non-linear evo- 
lution. This model is able to reproduce the matter power 
spectrum of V-body simulations over a wide range of scales 
and relates the linear and non-linear power at the same 
scale through fitting formulae. In a more recent development, 
Smith et al. (2002) have presented a new set of fitting func- 
tions based on the Halo Model and calibrated to a set of 
iV-body simulations. 

Barber & Taylor (2003) have now shown excellent agree- 
ment between the power spectrum in the lensing convergence 
obtained from numerical simulations in which they calcu- 
lated the full 3-dimensional shear matrices along lines of 
sight and the predictions from the Halo model fitting func- 
tions of Smith et al. (2002). 

An alternative approach used to describe the probabil- 
ity distribution function of the density field, from which the 
associated weak lensing statistics can be obtained, has been 
based directly on the many-body correlations. The most 
common model of this kind expresses the p— point corre- 
lations as a sum of p — 1 products over the two-point corre- 
lations linking all p points. This yields the class of "Hierar- 
chical models," which are specified by the weights given to 
any such topology associated with the p — 1 products (e.g., 
Fry 1984, Schaeffer 1984, Bernardeau & Schaeffer 1992, and 
Szapudi & Szalay 1993, 1997, Munshi et al. 1999a, Mun- 
shi, Melott, Coles 1999d). Once these weights have been as- 
signed it is possible to resum all many-body correlations 
and to compute the pdf of the density field, or of any quan- 
tity which is linearly dependent on the matter density. This 
is most easily done for "minimal tree-models", where the 
weight associated to a given tree-topology is set by its ver- 
tices (e.g., Bernardeau & Schaeffer 1992, Munshi, Coles & 
Melott 1999b, Munshi, Coles & Melott 1999d), or for "stel- 
lar models", which only contain stellar diagrams (Valageas, 
Barber & Munshi, 2003). 

Using such an approach, coupled to the Hamilton et al. 
(1991) prescription for the two-point correlation, Valageas 
(2000a, b) and Munshi & Jain (2000 and 2001), Munshi 
& Coles (2000,2002) were able to compute the pdf of the 
weak-lensing convergence whilst the associated bias was con- 
sidered by Munshi (2000) and the cumulant correlators as- 
sociated with such distributions were evaluated by Munshi 
& Jain (2000). Munshi & Wang (2003) further extended 
these studies to cosmological scenarios including dark en- 
ergy. These methods can also handle more intricate quanti- 
ties like the aperture-mass or the shear which involve com- 
pensated filters and require a detailed model for the many- 
body correlations. Thus, using a minimal tree-model for the 



non-linear regime, Bernardeau and Valageas (2000) were 
able to predict the pdf of the aperture-mass and to ob- 
tain a good agreement with numerical simulations (they 
also showed that similar techniques could be applied to the 
quasi-linear regime where the calculations can actually be 
made rigorous). On the other hand, adopting a stellar model 
for the many-body correlations, Valageas, Barber & Munshi 
(2004) obtained excellent agreement for the shear pdf when 
compared with the results of V-body simulations. In a more 
recent paper, Barber, Munshi & Valageas (2004) described 
the results for the pdf of the convergence based on hierarchi- 
cal models and again showed excellent agreement with the 
results from iV-body simulations. 

In this paper, we employ the minimal-tree model and 
the stellar model for the density field to predict statistics 
for the weak lensing shear and aperture mass, which cor- 
respond to a top-hat smoothing filter and a compensated 
filter, respectively. While earlier studies (e.g., Bernardeau & 
Valageas 2000; Valageas, Barber & Munshi 2004) showed 
that such approaches provide a promising tool to obtain 
quantitative predictions for weak-lensing observables, we 
compare in details in this article these theoretical predic- 
tions against results from numerical simulations, over a large 
range of scales and redshifts. Thus, we are able to check that 
these methods provide indeed reliable means to predict weak 
lensing statistics, from quasi-linear scales up to highly non- 
linear scales. In addition, we compare the predictions ob- 
tained from the minimal-tree model and the stellar model. 
This allows us to investigate the sensitivity of these weak- 
lensing observables onto the detailed angular behaviour of 
the many-body density correlations (their overall amplitude 
at a given scale being identical for both models). We find 
that such a dependence only shows up in the negative tail 
of the aperture- mass at non-linear scales (which can then be 
used to discriminate between different angular models while 
other regimes set the amplitude of the density correlations). 
Therefore, this paper extends to the weak lensing shear and 
aperture mass the detailed analysis presented in Barber et 
al.(2004) for the convergence (which is somewhat simpler). 

In section 2, we describe weak lensing distortions in gen- 
eral and the use of filters to smooth the data. In section 3, 
we introduce the minimal-tree and stellar models for the 
density field and in section 4, we describe how the pdfs for 
the weak lensing observables are derived in these models. 
In section 5, we describe the V-body simulations and the 
procedure for computing the lensing statistics in the simu- 
lations; we also outline the procedure for binning the data 
and applying the smoothing filters. In section 6, we declare 
the results in detail which are discussed finally in section 7. 



2 WEAK LENSING DISTORTIONS 

As a photon travels from a distant source towards the ob- 
server, its trajectory is deflected by density fluctuations close 
to the line-of-sight. This leads to an apparent displacement 
of the source and to a distortion of the image (as the de- 
flection varies with the direction on the sky). These effects 
can be observed through the amplification or the shear of 
the observed images of distant sources. One such measure 
of the distortion may be given in terms of the convergence 
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along the line-of-sight, given by (e.g., Bernardeau et 

al., 1997; Kaiser, 1998): 



k{6) 



with: 



3^ n 



dx w(x,Xs) <5(X>W)> 



c 2 ©(Xs) 



X) 



(1 + *). 



(1) 



(2) 



where z corresponds to the radial distance x an d 2? is the 
angular distance. Here and in the following we use the Born 
approximation which is well-suited to weak-lensing studies: 
the fluctuations of the gravitational potential are computed 
along the unperturbed trajectory of the photon (Kaiser, 
1992). Thus the convergence k(#) is merely the projection 
of the local density contrast 8 along the line of sight up to 
the redshift z a of the source. Therefore, weak lensing obser- 
vations allow us to measure the projected density field k(#) 
on the sky. Note that by looking at sources located at differ- 
ent redshifts one may also probe the radial direction. From 
eq.(l) we can see that there is a minimum value, K m i n (z s ), 
for the convergence of a source located at redshift z s , which 
corresponds to an "empty" beam between the source and 
the observer (8 = — 1 everywhere along the line of sight): 



30n r 

2 Jo 



dx w(x,Xs)- 



(3) 



Following Valageas (2000a, b) it is convenient to define the 
"normalized" convergence, k, by: 



|Kmin| J 



dx w 8, with w 



J Xs d X w{x,Xs) 



■ (4) 



which obeys k > — 1. Here we introduced the "normalized 
selection function," w(x,Xa)- 

In practice, one does not study the projected density 
k(i?) itself but applies first a smoothing procedure. For in- 
stance, it is customary to investigate the "smoothed conver- 
gence" k s (where the subscript "s" refers to "smoothed") 
given by: 



k s =J 



with: 



8Q9 < 9 S 



(5) 



(6) 



where 0($ < 9 S ) is a top-hat with obvious notations. Thus, 
the "smoothed convergence" k 3 is simply the projected den- 
sity contrast k($) smoothed with a normalised top-hat U K 
of angular radius 8 3 . By varying the radius 6 S one probes 
the density field at different wavelengths. 

The "smoothed convergence" k s was already studied at 
length in previous papers (e.g., Barber, Munshi & Valageas 
2004). Here we focus on the shear 7 and the aperture mass 
M a p. These two quantities can again be expressed in terms 
of the projected density contrast k(i?) as in eq.(5). Thus, the 
smoothed shear 7 S = 7i s + iy^s corresponds to the filter XJ 1 
given by (see Valageas, Barber & Munshi 2003): 



t/ 7 (0) = - 



e(fl > e a ) l2a 

7T0 2 6 ' 



(7) 



where 0(i? > 6 S ) is again a Heaviside function with obvious 
notations and a is the polar angle of the vector 1? with the 



1-axis (the 2-axis has a = 7r/2). The comparison of eq.(7) 
with eq.(6) shows that while the smoothed convergence only 
depends on the matter within the cone formed by the angular 
window, 6 S , the smoothed shear only depends on the matter 
outside this cone. Then, the smoothed shear component 71 
along the 1-axis is described by the filter U 11 : 



u * 9(0 > 9s) 

U~!lW = ^2 COS2Q - 



(8) 



Note that the filters U y and U 11 depend on both the length 
and the angle of the two-dimensional vector $ in the plane 
perpendicular to the mean line-of-sight (in this article we 
only consider small angular scales). Moreover, they are com- 
pensated filters which is an interesting feature since conver- 
gence maps are only reconstructed up to a mass sheet de- 
generacy. In fact, the shear 7 is actually the quantity which 
is most directly linked to observations. 

One drawback of the shear components is that they are 
even quantities (their sign can be changed through a rota- 
tion of axis, see eq.(7)) hence their third-order moment van- 
ishes by symmetry and one must measure the fourth-order 
moment (7^ (i.e. the kurtosis) in order to probe the de- 
viations from Gaussianity. Another quantity which involves 
a compensated filter but is not even is the aperture-mass 
M a p. It simply corresponds to a compensated filter with po- 
lar symmetry: t/M ap (#) = Vm„ v ($)• For instance, following 
Schneider (1996) one can use: 



6(i? < 



9 



1-^ 

fl2 



d 2 



(9) 



The advantage of such compensated filters is that one can 
also express M ap as a function of the tangential component 
7t of the shear (Kaiser et al. 1994; Schneider 1996) so that 
it is not necessary to build a full convergence map from 
observations. Besides, the aperture-mass provides a useful 
separation between E and B modes. 

In the following, we shall write U$ for arbitrary filters 
when our results apply to any quantity X defined as in eq.(5) 
with some filter U%, and we shall specify U K , ..,Um &p for 
particular cases. For some purposes it is convenient to work 
in Fourier space, therefore we define the Fourier transform 
of the density contrast by: 



*(x) 



dk e <k x <5(k), 



(10) 



where x and k are comoving coordinates. Then, we define 
the power-spectrum P(k) of the density contrast by: 



(<5(ki)<5(k 2 )> = Mki + k 2 ) P(ki), 



(11) 



where 8d is Dirac's distribution. This yields for the two- 
point correlation £2(2:) of the density contrast: 

£ a (a0 = (J(xi)J(xi + x)> = J dk e lk ' x P(k). (12) 

We also introduce the power per logarithmic interval A 2 (fc) 



A 2 (k,z) =4Tvk 3 P(k,z). (13) 

Then, we can write for any normalised quantity X (such as 
the normalised smoothed shear -%): 



X 



-J- 
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J d X w J < 



dx w / dk e lfc ii x (kx2?6> s ) <5(k), (15) 



where fey is the component of k parallel to the line-of-sight, 
kx is the two-dimensional vector formed by the components 
of k perpendicular to the line-of-sight and W% is the Fourier 
form of the real-space filter U ^ : 



(16) 



Following previous works we explicitly introduced the angu- 
lar scale 9s in the definition of the Fourier filter . For the 
smoothed convergence k s we obtain: 

2j 1 (k ± ve s ) 



W K Qt±D0.) = 



k ± V8 s 



(17) 



where Ji is the Bessel function of the first kind of order 1, 
while for the shear we have: 



W~ f (k ± V6 s ) = W K (k±V0 3 ) e l2a 
and: 

W 71 (kxD<9 s ) = W K {k±T>0„) cos 2a, 



(18) 



(19) 



where a is the polar angle of the transverse wavenumber 
kx- Finally, we obtain for the aperture-mass defined from 
the filter (9): 



2AJ i {k ± Ve s 
(k ± V6s) 2 



(20) 



In any case, the variance of the observable X is obtained 
from eq.(15) as: 



{X\ = / d X w 

10 



f . a 1 f dk_x A 2 (k 

2 



,, i ., k± w^ve.y. (2i) 



3 THE DENSITY FIELD 

In order to derive the properties of weak lensing observables, 
like the shear 7, we clearly need to specify the properties 
of the underlying density field. In this paper we use sim- 
ple models which we described in detail in Valageas, Barber 
& Munshi (2003) and Barber, Munshi & Valageas (2003), 
hence we only briefly recall the main elements of this frame- 
work here. First, we describe the pdf V(Sr) of the density 
contrast at scale R through its generating function tp(y): 



/DC 
-1 



d8 R e 



(22) 



where Sr is the density contrast within spherical cells of 
radius R and volume V while £ 2 is its variance: 



)>: j ^ <5(x) and £ 2 = (6%). 



(23) 



The function tp(y) defined from eq.(22) is also the generating 
function of the cumulants of the density contrast and its 
expansion at y — reads: 



<p(y) = 



S p y p with S p = 1 



p-1 



(24) 



As described in Valageas, Barber & Munshi (2003) and Bar- 
ber, Munshi & Valageas (2003) we parameterize the gener- 
ating function ip(y) through the skewness S3 of the density 
field, which we estimate as: 



(25) 



Here, S^ L is the exact result derived in the quasi-linear limit 
while is the prediction of HEPT (Scoccimarro & Frie- 
man 1999) for the highly non-linear regime. We also intro- 
duced the density contrast A v i r at virialization which marks 
the onset of the highly non-linear regime. Eq.(25) applies to 
intermediate scales. At large scales (A 2 (k s ,z) < 1) we sim- 
ply take 53 = S® h while at small scales (A 2 (k s ,z) > A vir ) 
we use S3 = Ss L . Finally, we define the typical wavenumber 
k s {z) probed by weak lensing observations as: 



1 4 

ks( z ) = for k or 7; k 3 (z) — 



V(z)6 s 



V{z)6 a 



for M ap . (26) 



Indeed, as seen from Fig. 3, for the same angular radius 
9 a the aperture-mass M ap probes higher wavenumbers than 
the convergence or the shear. This is obviously due to the 
radial structure of the filter [/j( ap (t?). As seen from eq.(22), 
in order to determine the pdf V(5r) we simply need to add a 
prescription for the variance £ 2 . We use the fit to numerical 
simulations provided by Peacock & Dodds (1996). 

^From the pdf P(8r), or the generating function ip(y) 
and the variance £ 2 , we have a full description of the den- 
sity fluctuations smoothed over spherical cells. As shown in 
detail in Valageas (2000) and Barber, Munshi & Valageas 
(2003) this is sufficient to obtain the properties of the 
smoothed convergence k s . However, as noticed in Bernardeau 
& Valageas (2000) and Valageas et al. (2003), this informa- 
tion is no longer sufficient when we study more intricate ob- 
servables like the shear or the aperture-mass which involve 
compensated filters. Therefore, we need to specify the de- 
tailed angular behaviour of the many-body correlation func- 
tions £ p (xi, .., Xj,), defined by (Peebles 1980): 



£p(xi,..,x p ) = (<5(xi)..<$(x p )) c 



(27) 



As in Barber, Munshi & Valageas (2003) we shall compare 
two specific cases within the more general class of "tree- 
models" . The latter are defined by the hierarchical property 
(Schaeffer, 1984, and Groth & Peebles, 1977): 



(28) 



(a) 



p-1 



p=2 



where (a) is a particular tree-topology connecting the p 
points without making any loop, Q p a ^ is a parameter as- 
sociated with the order of the correlations and the topology 
involved, t a is a particular labeling of the topology, (a), and 
the product is made over the (p — 1) links between the p 
points with two-body correlation functions. 

Then, the "minimal tree-model" corresponds to the spe- 
cific case where the weights Q p a ^ are given by (Bernardeau 
& Schaeffer, 1992): 

Q P a) = n v " ( 29 > 

vertices of (a) 

where v q is a constant weight associated to a vertex of the 
tree topology with q outgoing lines. The advantage of this 



© 0000 RAS, MNRAS 000, 000-000 



Statistics of weak lensing shear and aperture-mass 5 



minimal tree-model is that it is well-suited to the compu- 
tation of the cumulant generating functions as defined in 
eq.(24) for the density contrast Sr. Indeed, for an arbitrary 
real-space filter, F(x), which defines the random variable s 
as: 



/ 



dxF(x)5(x) and £ s = (s 2 ), 



(30) 



it is possible to obtain a simple implicit expression for the 
generating function, ip s (y) (see Bernardeau & Schaeffer, 
1992, and Jannink & Des Cloiseaux, 1987): 



<Ps(y) 



r(x) 



dx F(x) 



C[r(x)]- 



-WC[r(x)l 



/|dx' F(x') £[r(x')] 



(31) 



(32) 



where the function Cf( r ) is defined as the generating func- 
tion for the coefficients v v : 



p=i 



(-If 
p! 



with 



(33) 



A second simple model is the "stellar model" introduced 
in Valageas, Barber & Munshi (2003) where we only keep 
the stellar diagrams in eq.(28). Thus, the p— point connected 
correlation £ p of the density field can now be written as: 

p 

The advantage of the stellar-model (34) is that it leads to 
very simple calculations in Fourier space. Indeed, eq.(34) 
reads in Fourier space: 

(5(k 1 )..5(k p )) c =^fo(k 1 + ..+k p )^J|P(fc,). (35) 

i — l j^i 

Following Valageas, Barber & Munshi (2003) and Barber, 
Munshi & Valageas (2003) we shall use the simple approxi- 
mation S p ~ S p . Alternatively, we may define the function 
<p(y) as the generating function of the coefficients S p , rather 
than S p , through its Taylor expansion at y = 0. 



4 THE PDF OF WEAK LENSING 
OBSERVABLES 

^From the models described in Sect. 3 for the density field 
we can derive the pdf of the weak lensing observables defined 
in Sect. 2. The procedure is described in detail in Barber, 
Munshi & Valageas (2003) (and references therein). Hence 
we only briefly recall here the main results. 

4.1 Minimal tree-model 

We first apply the minimal tree-model to an arbitrary weak- 
lensing observable X defined by the filter . As described 
in Valageas (2000b), in order to derive the pdf V(X) we first 
compute the cumulants {X p ) c , which we resum to obtain the 
generating function fxiv) defined as in eq.(24): 



'fix 



(»)=E 



P \ (xYc' 1 



(36) 



Next, the pdf V(X) is given by the inverse Laplace trans- 
form: 



/ + IC 
=, e 
. ioo 27Ti(V2) c 



Ay r [Xy- Vi ,(y)]/(X 2 ) c _ ^ 

Thus, we first obtain the cumulants (X p ) c from (14): 



J° i=2 * i=l 



x Cp ( ' X? 2 '-' v$p ;z ) ■ (38) 

Then, as seen in Valageas (2000b) we note that within the 
framework of a tree-model (28) the 2-d correlations u> p ex- 
hibit the same tree-structure, with: 



-OO ■ n \ 



Xp 

v3 v 



(39) 

In terms of the power-spectrum P(k) we can also write ui2 
as: 

dfc A 2 {k, z) 
k k 



W2(l?l, $2 5 z) = TV 



Jo{KD\di-$2\). (40) 



Then, in the case of a minimal tree-model (29) we can per- 
form the resummation (31)-(32) for the 2-d correlations ui p , 
since the latter obey the same minimal tree-model. This 
yields (see Bernardeau & Valageas 2000 and Barber, Munshi 
& Valageas 2003 for details): 



Vcyi. yw 



(41) 



10 -a* V <V 2 ) C 

where we introduced the 2-d generating function tp cy i. associ- 
ated with the 2-d correlations u> p , given by the resummation: 



VcyUy) = y d0Ux(0) 



(42) 



-y J dS'U^') "i^\ Z) C[r(<)')] ( 13 



Here we introduced the angular average lo 2 -^ of the 2-d cor- 
relation ll>2 , associated with the filter {7 X : 



/ 



"2x(z) = / dtfidtfa U^i)U^2) D2(&i,tir,z). 



(44) 



Thus, we obtain in this way both generating functions tpj le 
and ip M ~ associated with the smoothed normalised shear 

component 7i s and the aperture- mass M ap . This yields in 
turn the pdfs Vtfu) and P(M ap ), using eq.(37). Note that 
for the aperture-mass the implicit system (42)-(43) simpli- 
fies somewhat since the filter Um^ ($) only depends on the 
length 

4.2 Stellar model 

Next, we can use the same procedure within the framework 
of the stellar model (34). Thus, working in Fourier space we 
obtain from eq.(15) the cumulants (X p ) c as: 

(x p )c = f ^(2™) P / ffcikx, w A (k ±j ve s ) 
J J j=i 
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x S p fo(k ± i + .. + k ±p )P(k ±2 )..P(k ±p ). (45) 

Next, using the standard exponential representation of the 
Dirac distribution (see Valageas, Barber & Munshi 2003) 
and eq.(16), we can write: 

(X p ) c = jT' d X w p J M U A 0) S P IAxJT' 1 , (46) 



where we introduced: 

Then, using eq.(24) we obtain: 



Wx(k_LT>e s )e ik±/D9 . (47) 



<Px(y) 



= T dx I 



<p yw- 



■z\. (48) 



The result (48) allows us to obtain both generating functions 
tpj ls and ifi M ~ associated with the smoothed shear compo- 
nent 7i s and the aperture-mass M ap . This again yields in 
turn the pdfs P(7i s ) and P(M ap ), using eq.(37). For the 
aperture- mass eq.(48) also simplifies somewhat since the fil- 
ter L?M ap (#) and Im^ p (x^) only depend on the length 

4.3 Exponential tails 

As described in Bernardeau & Schaeffer (1992), the implicit 
system (42)- (43) usually yields branch cuts along the real 
axis for the generating function f C y\.(y)- In fact, we actually 
define the generating function <p(y) of the 3-d density con- 
trast through a similar implicit system (see Valageas, Barber 
& Munshi 2003 and Barber, Munshi & Valageas 2003 for de- 
tails), so that the function <p(y) shows a branch cut along 
the negative real axis for y < y s , with: 

«" = U+2J Wlth (49) 

Here S3 is the skewness of the density contrast at the scale 
and time of interest. The singularity y s leads to an expo- 
nential tail for the pdf P(5r) for large positive Sr. On the 
other hand, for large y the generating function ip(y) shows a 
slow power-law growth (see Bernardeau & Schaeffer 1992): 
1- 



y — > +00 : ip(y) + y ~ y w with < u> < 1. 



(50) 



This large y behaviour leads to a strong cutoff at low den- 
sities and V(5r) = for 8r < —1 (this lower bound is 
set by the coefficient of the term linear over y in eq.(50)) 
when we can push the integration path to +00 in eq.(37). 
As seen in Valageas (2000) and Barber, Munshi & Valageas 
(2003), since the smoothed convergence is described by a 
2-d top-hat, which is quite similar to the 3-d top-hat associ- 
ated with the density contrast 8r, the generating functions 
<Pcy\.{y), <Pk B and the pdf V(k s ) show the same behaviour 
as for Sr. In fact, as shown in Valageas (2000) and Barber, 
Munshi & Valageas (2003), one can directly use ip(y) and 
V(5r) to obtain up to a good accuracy the properties of the 
smoothed convergence. 

Of course, for more intricate filters the properties of 
the generating function f%(y) and of the pdf V{X) can ex- 
hibit very different behaviours. In particular, for the shear 
component 7i s the generating function ip lls (y) and the pdf 
P(7i s ) are now even. This property is obviously preserved 
within both models used in this paper. For the minimal 



tree-model, this appears in the implicit system (42)-(43) 
through the factor cos 2a of the filter U yi (1?) given in eq.(8), 
which clearly implies that the function <p cy i. is even. On the 
other hand, for the stellar model this property shows up in 
cq.(48) because the angular dependence of the filter U lt {•&) 
leads to 7 71 (x, oc cos 2a, see Valageas, Barber & Munshi 
(2003) for details. Then, in both cases the generating func- 
tion <fi lls (y) shows two symmetric branch cuts along the real 
axis for y < -y 3 , lle and y > y s , lls (with y s>lls > 0). For 
the stellar model, this singularity can be expressed in terms 
of y s as (see Valageas et al. 2003): 



y s , ll3 = mm 



th- 



ills). 



(51) 



Then, the pdf Vtfu) shows two symmetric exponential tails 
for jis — > ±00. 

On the other hand, the aperture-mass M ap can lead to 
a more intricate behaviour. Since it involves a compensated 
filter like the shear component 7i s , we can expect extended 
tails both for large positive and negative M ap . However, the 
generating function ip M ~ (y) and the pdf P(M a . p ) are not 
even. From the shape of the filter (9) we can actually expect 
the fall-off in the pdf to be sharper for negative M ap . For 
the stellar model, the two branch cuts along the real axis of 
the generating function ip M - (y) are given by: 



y, 



JJs 



' Map max z ^(w/Af ap /(M ap 2 ) c ) 



y~ M ~ < 0, (52) 



and 

-4- 

y 



9 — , y\ f > 0. (53) 



min Zi #(w/M ap /(M ap ) c ) 



< V .,- which 



We can check numerically that \y 

means that the cutoff of V(M !!l p) is indeed stronger for neg- 
ative M ap . This property is also verified by the tree-model. 

Here, we must point out that such exponential tails for 
the various pdf (associated with a branch cut along the neg- 
ative real axis for ip(y)) are a mere consequence of the pa- 
rameterization we use for (p(y). Indeed, as recalled above, 
the pdf V(k s ) of the smoothed convergence closely follows 
the pdf V{5r) of the 3-d density contrast. Therefore, if we 
model V(Sr) in such a way that it shows a stronger cutoff 
than the exponential at large overdensities (such as a Gaus- 
sian) it would translate into a similar behaviour for V(n s ) 
and the generating functions <p{y) and <p K (y) would have no 
branch cut. Such a behaviour would also apply to the shear 
or the aperture-mass. Similarly, a cutoff which is smoother 
than exponential for V(5r) would yield extended tails for 
weak lensing observables. In the quasi-linear limit it is pos- 
sible to derive the tails of the pdf T(5r) and to show that the 
high-density cutoff actually involves the exponential of some 
power-law (see Valageas 2002a) but this discrepancy is negli- 
gible in the range of interest (the far tail corresponds to very 
rare events which are irrelevant for most practical purposes). 
In the non-linear regime there is no rigorous derivation of the 
tails of the pdf V(Sr) (however see Valageas 2002b for more 
details) but the simple phenomenological model recalled in 
Sect. 3 appears to be sufficient, as seen in Barber, Munshi 
& Valageas (2003) from a comparison with numerical simu- 
lations for the weak lensing convergence n s . 
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4.4 The tail of"P(M ap ) for negative M ap at small 
angles: sensitivity onto the angular behaviour 
of the many-body correlations 

We can note that at small scales, 8 S <, 0.4', the factor 
7~M ap becomes positive for all z and which means that 
the branch cut for positive y is repelled to +00 within the 
stellar model. Then, the pdf V(M a . p ) actually vanishes for 
Map < 0. Since we always have (Map) = 0, this means that 
the pdf 7 , (M ap ) actually contains a singular part at the ori- 
gin. This behaviour appears when the power per logarithmic 
interval A 2 (k) is natter than k. For CDM-like power-spectra 
this corresponds to small scales (i.e. to 6 a Ss 0.4', see Fig. 3). 
Note that this behaviour is independent of the coefficients 
S p , that is of our parameterization for the generating func- 
tion <p(y) or the pdf V(Sr). It is only due to the angular 
structure implied by the stellar model and would remain un- 
changed whatever the prescription used for the pdf V(5r) of 
the density contrast. This suggests that no physical density 
field can be exactly described by the stellar model (34). On 
the other hand, the minimal tree-model does not show such 
a breakdown at small angles and it always yields a smooth 
exponential tail for large negative M ap . 

This discrepancy between both models shows that the 
tail of the pdf 7 ? (M ap ) for negative M ap is very sensitive to 
the detailed angular properties of the density field. This is 
an interesting feature since it means that one could extract 
in principles some useful information about the structure of 
the density field at small scales from the aperture- mass M ap . 
Unfortunately, such small angular scales may be difficult to 
study in actual observations. However, as seen in Sect. 6.6 a 
large discrepancy between both models is already apparent 
for a smoothing radius 8 S — 2' which may be within the 
reach of observations. 



4.5 Integration path 

In order to compute the pdf V(X) from the inverse Laplace 
transform (37) we must perform the integration over y in 
the complex plane. We choose the integration path so that 
the argument of the exponential in eq.(37) is a real negative 
number in order to avoid oscillations and to obtain a fast 
convergence. From the definition (36) it is clear that the 
generating function <p%(y) obeys the symmetry <Px(v*) = 
< Px{v)*j f° r an y rea l observable X like the shear component 
7i s or the aperture-mass M ap , so that the integration path 
over y is symmetric with respect to the real axis. Moreover, 
we have to make sure that the integration path does not 
cross the branch cuts of ip^{y). This implies that for large 
\X\ the integration path is pinched on the real axis onto the 
singularities y This directly yields an exponential tail 

for the pdf V(X), as can be seen in a straightforward way 
from the expression (37). We discussed this point in more 
detail in Sect. 4.3 above. For illustration, we show in Fig. 1 
the integration paths obtained for the aperture-mass in the 
LCDM case for z s = 1 and 9 S = 4' '. For M ap = the path 
runs through the origin y — while for large positive M ap 
it gets stucked onto the negative branch cut of tp^ (y) (and 
similarly for large negative M ap ). 



0.2 



>: 



M = o\ 

ap \ 


1 

M = -0.01/ 




z,= l, 0„=4' 

1 



-0.2 





y r 



0.2 



Figure 1. The integration path over the complex plane for the 
aperture-mass M ap for the LCDM scenario, the source redshift 
z s = 1 and the angular radius 8 a =4'. We show the paths associ- 
ated with M ap = -0.01 (right), A/ ap = (center) and M ap = 0.01 
(left). We clearly see that for large |M ap | the integration path gets 
stucked onto one of the branch cuts of the generating function 

vm- GO- 



4.6 Edgeworth expansion 

As is well known, a centered Gaussian distribution is fully 
defined by its second order moment {X 2 ) c . When the pdf 
V(X) deviates from the Gaussian all higher-order cumulants 
{X p ) c generically become non-zero and the pdf now depends 
on this whole series. Thus, in Sect. 4.1 and Sect. 4.2 we had 
to resum all these cumulants through the generating func- 
tion fxiv) d enne d in eq.(36) in order to determine the pdf 
V{X) as in eq.(37). However, when the deviations from the 
Gaussian are small one can use the asymptotic Edgeworth 
expansion which describes the departures from the Gaussian 
encoded by the few lowest order cumulants. Let us define the 
parameters Sp associated with the variable X as in eq.(24): 



p=2 



with S„ 



{X p )c 



2\P-i 



(54) 



Here we assumed the random variable X to have zero mean 
({X) = 0) as is the case for weak lensing observables from 
eq.(14) since (5) = 0. Then, substituting the expansion (54) 
into eq.(37) and expanding the non-Gaussian part of the 
exponent one obtains (e.g., Bernardeau & Kofman 1995): 



V{X) 



si 

24 



2 /2 



HAv) + 



72 



-H 6 (u) 



with: 

o x = { xyj 2 



and v — 



X 



(55) 



(56) 
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Here we introduced the Hermite polynomials H n (v). In par- 
ticular we have: 

Hz(v) = v 3 - 3u and H A (y) = v A - 6u 2 + 3. (57) 

Since this is an asymptotic expansion, the Edgeworth ex- 
pansion (55) is only useful for moderate deviations from the 
Gaussian, that is when the first correcting term is smaller 
than unity (typically \u\ £ 1 and |cr J fS'^ f | ^ 1), and the ac- 
curacy does not improve by including higher order terms in 
the expansion (55) (unless the pdf is extremely close to the 
Gaussian). The advantage of the Edgeworth expansion is 
that it provides a straightforward estimate of the pdf V(X) 
when it is still relatively close to the Gaussian, using only 
the skewness S$ or the kurtosis S 4 . This avoids the need 
to resum all higher-order cumulants and to perform the in- 
tegration over the complex plane as in eq.(37). 



5 THE SHEAR AND APERTURE-MASS 
STATISTICS FROM NUMERICAL 
SIMULATIONS 

The numerical method for the computation of the lensing 
statistics is based on the original formalism of Couchman, 
Barber &: Thomas, 1999, and which has been further de- 
veloped by Barber, 2002. The original formalism allows for 
the computation of the three-dimensional shear matrices at 
a large number of locations within each iV-body simulation 
volume output. In the present work we have evaluated the 
shear at 300 locations along every one of 455 x 455 lines of 
sight in each of the simulation volumes. The development al- 
lows for the successive combination of these matrices along 
the lines of sight and throughout the linked simulation vol- 
umes from the sources at the required redshifts to the ob- 
server at z — 0. The overall procedure therefore gives rise to 
Jacobian matrices for each line of sight and for each of the 
specified source redshifts from which the required lensing 
statistics have been generated. 

Our procedure has been applied to two different cos- 
mological simulations created by the Hydra Consortium* 
using the 'Hydra' iV-body hydrodynamics code (Couchman, 
Thomas & Pearce, 1995). The parameters describing the two 
cosmologies, LCDM and OCDM, are given in Table 1. Both 
contained 86 3 dark matter particles of mass 1.29 x lO 11 ^ 1 
solar masses each, where h is the value of the Hubble pa- 
rameter expressed in units of 100 km s _1 Mpc -1 . We used 
a variable particle softening within the code, to reflect the 
density environment of each particle, the minimum value of 
which, for particles in the densest environments, was chosen 
to be 0.0007(1 + z) in box units, where z is the redshift of 
the particular simulation volume. 

The angular resolution in the LCDM cosmology was 
0'.34, equivalent to the minimum value of the particle soft- 
ening at z — 0.36, which is the redshift for the maximum 
lensing effects for sources at a redshift of 1 in that cosmol- 
ogy. In the case of the OCDM cosmology, the angular reso- 
lution was 0'.37. The angular size of the survey referred to 
in Table 1 corresponds to completely filling the front face of 
the redshift 1 simulation volume for an observer at redshift 

* (http://hydra.mcmaster.ca/hydra/index.html) 



Table 1. The parameters used in the two cosmological simula- 
tions. f2 m is the matter density parameter, Qj^ is the vacuum 
energy density parameter, T is the power spectrum shape param- 
eter, ag is the normalisation on scales of 8/i _1 Mpc, 9 rcs is the 
angular resolution and d^ urvcy is the angular size of the complete 
field of view throughout the simulations. 





S7 m 


n A 


r 


0"8 


$rcs 


1 


? 2 

'survey 


LCDM 
OCDM 


0.3 
0.3 


0.7 
0.0 


0.25 
0.25 


1.22 
1.06 


0'.34 
0'.37 


2° 
2° 


.6 x 2°. 6 
.8 x 2°. 8 



zero. For source redshifts greater than 1, the periodicity of 
the particle distributions was used to allow lines of sight be- 
yond the confines of the simulation volumes to be included, 
as described in Barber, Munshi & Valageas, 2003. 

The sources for our simulations were assumed to lie in 
a regular 455 x 455 array at the front face of the simulation 
volume corresponding to the specified redshift. We used 14 
different values for the source redshifts in each cosmology, 
whose exact redshift values are specified in Table 2. The 
redshifts selected were chosen to be close to redshifts of 0.1, 
0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 3.0 and 3.5. 
In this paper the redshifts are referred to loosely as these 
approximate values, although in the determination of the 
lensing statistics, the actual redshift values were used. 

Each of the simulation volumes had comoving side- 
dimensions of 100fe _1 Mpc and to avoid obvious structure 
correlations, each was arbitrarily translated, rotated and re- 
flected about each coordinate axis for each of the total of 
N = 10 runs through all the volumes. In this way, we ob- 
tained N sets of Jacobian matrices for the different runs. 

By extracting the data needed from the Jacobians, 
the data for the shear components and convergence were 
smoothed on the different angular scales using a top-hat fil- 
ter for the shear and the compensated filter for the aperture 
mass. To obtain the required pdfs and lower order moments 
the unsmoothed field in real space was convolved with the 
appropriate filters using the following scheme. Concentric 
annuli, separated equally in the radial direction, were used 
for sampling of the denisty field, as shown in figure 2. These 
annuli were than divided into equally spaced angular bins. 
The density field was computed along the radial and angular 
bins using a linear interpolation scheme from adjacent bins, 
before convolving it with the filters. Various combinations 
of the radial bin width and angular bin size were consid- 
ered to check the convergence and numerical stability of the 
scheme. Our scheme enforces the radial symmetry inherent 
in the window functions which is important when dealing 
with window functions that have more features compared 
to a tophat filter. We have checked the stability by taking 
up to 20 radial binning divisions for large smoothing radius 
and 60 divisons along the angular directions. In addition, 
we have given random shifts to the rectangular grid along 
both axes to get a better sampling of the density field. We 
have also checked various levels of dilution by considering 
randomly selected points in the catalogue. In addition, we 
have checked various levels of uniform dilution where we 
choose uniformly spaced points to evaluate the statistics of 
the aperture mass, M ap . It became clear from our studies 
that grid effects will be much more pronounced in any sta- 
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Table 2. The redshifts of the sources in the two cosmologies. 
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Figure 2. Schematic description of the angular bins used for 
evaluation of M ap statistics. 

tistical analysis related to compensated niters. This is due 
to the fact that for a given smoothing scale, compensated 
filters take more contributions from smaller angular scales. 

Finally, the computed values for the pdfs, higher-order 
moments and aperture mass from each of the runs in each 
cosmology were averaged so that the errors on the means of 
la/yN for each statistic were determined. 



6 RESULTS 

Our study can be divided in two parts. In addition to the 
lower order moments we compute the complete pdf and we 
study their variation as a function of smoothing radius as 
well as their evolution with redshift. For the pdf of the aper- 
ture mass we investigate both the minimal tree model and 
the stellar model, which allows us to study the dependence 
on the detailed modeling of the correlation function (both 
models give close results for the shear). 

6.1 The length scales probed by shear statistics 

As we study a range of source redshifts and smoothing an- 
gular scales the weak lensing observables (such as the shear 
or the aperture-mass) probe the density field over various 
length scales and redshifts, which run from the linear to 
the highly non-linear regime. The redshift dependence of 
the typical length scale comes from the angular diameter 
distance T>(z), see eq.(26). The typical comoving wavenum- 
ber probed by various shear statistics for a smoothing angle 
6s = 1' are of the order of 1 to lOftMpc -1 which corre- 
sponds to scales 0.1 to l/i" 1 Mpc, see Fig. 5 in Barber, Mun- 
shi & Valageas (2003). In terms of the power per logarith- 
mic wavenumber interval A 2 (k), for 1' < 6 S < 8', this yields 
0.1 <, A 2 (k) <: 400, see Fig.6 in Barber, Munshi & Valageas 




- 1 1 

log(k) [ Mpc-' ] 



Figure 3. The contribution of various comoving wavcnumbers k 
to weak lensing observables at redshift z = 1 along the line of 
sight. The solid line is the non-linear power A 2 (k)/k (obtained 
from Peacock & Dodds 1996) while the dotted line is the linear 
power A 2 L (k)/k (i.e. using the linear power-spectrum). The left 
dashed line shows the contribution W K (k±'D8 a )A 2 (k)/k to the 
variance of the smoothed convergence k s (or to the smoothed 
shear component 7i s ) for the angular radius d a = 1'. The right 
dashed line shows the contribution WM^ p {ki_Ti9 3 )A 2 (k)/k, mul- 
tiplied by a factor 5, to the variance of the aperture-mass M ap 
for 8 S = 1'. The upper axis shows the angular scale $ = l/T>k 
associated to comoving wavenumber k. 

(2003). Hence most of the contribution to the weak lens- 
ing shear comes from intermediate length scales which can- 
not be described by the quasi-linear limit nor by the stable- 
clustering ansatz. As a consequence, in order to study weak 
lensing over the angular scales of interest for observational 
purposes (i.e. from 1' up to 20') one needs to use a model 
which can be applied from the quasi-linear regime up to 
the highly non-linear regime. The models which we recalled 
in Sect. 3 (see also Valageas, Barber & Munshi 2003 and 
Barber, Munshi & Valageas 2003 for details) provide such 
a tool and yield all properties of the density field (both the 
amplitude and the angular dependence of the many-body 
correlations) over the entire dynamical range. 

We show in Fig. 3 the contribution of various comoving 
wavenumbers k to weak lensing observables at redshift z — 1 
along the line of sight. We first plot the non-linear power 
A 2 (k)/k (obtained from Peacock & Dodds 1996, solid line) 
and the linear power A 2 L {k)/k (i.e. using the linear power- 
spectrum, dotted line). Indeed, because of the projection 
associated with the integration along the line-of-sight, the 
power coming from a wavenumber k to weak lensing ob- 
servables is not A 2 (k) (as for the 3-d density contrast 5r) 
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Figure 4. Variance of the smoothed shear component ji s , 
(7? 3 ) 1 / 2 , as a function of the source redshift z 3 . The smooth- 
ing angle 8 3 is fixed at 1 arcminute in left panel and at 4 ar- 
cminute in the right panel. Lines correspond to the analytical 
prediction (21) and data points represent results from numerical 
simulations. Error bars are computed from scatter among various 
realizations. Variance of shear increases with source redshift z 3 
for a fixed smoothing angle. Dots represent LCDM simulations 
whereas OCDM simulations are represented by open circles. 
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Figure 5. Variance is plotted as a function of smoothing angle for 
various fixed redshifts. Left panel correspond to LCDM and the 
right panel correspond to OCDM cosmology. The solid (LCDM) 
and dashed (OCDM) lines in each panels denote the analytical 
prediction (21) for the variance at redshifts z 3 = 0.5, 1 and 1.5 
(from bottom to top). Points with error bars are numerical results 
from simulations. 



but A' 2 (k)/k, as seen from eq.(21). For CDM-like power- 
spectra, which flatten at small scales, we can see that the 
power A 2 (k)/k decreases beyond 6 Mpc — l . Next, we show 
the contribution W K (k±T>9 3 )A 2 (k)/k to the variance of the 
smoothed convergence k 3 (or to the smoothed shear compo- 
nent 7i s ) for the angular radius 9 3 = 1' (left dashed line). 
At small wavenumbers it follows the power A 2 (k)/k since 
W K (0) = 1 and it exhibits a cutoff beyond 1/(V9 3 ). The 
right dashed line is the contribution WM ap (k±V9 3 )A 2 (k)/k 
to the variance of the aperture-mass M ap for the same angu- 
lar radius 9 S = 1', multiplied by a factor 5 (for clarity on the 
figure). We can see that it probes a much narrower range of 
wavenumbers since the contribution from long wavelengths 
is suppressed because it involves a compensated filter with 
WM a p (0) = 0. Moreover, for the same angular radius we can 
check that M ap probes higher wavenumbers than k 3 because 
of this suppression of long wavelengths and of the profile of 
the filter which varies over scales ~ 9 3 /4. More precisely, 
from the upper axis which shows the angular scale $ = 1/T>k 
associated with comoving wavenumber k, we see that the 
smoothed convergence or the smoothed shear mainly probe 
the density field at the wavenumber k 3 = 1/T>9 3 while con- 
tributions to the aperture- mass peak around k 3 = 4/T>9 3 . 
This is why we chose the normalization (26) for the typical 
wavenumber k 3 . Finally, the curves in Fig. 3 clearly show 
that the weak lensing signal comes from non-linear scales 
which cannot be described by the linear power A\ (k)/k. 



6.2 Variance of Shear and Aperture-Mass 

We study the variance of smoothed shear components (7 2 s ) 1/ ' 2 
and aperture-mass (M ap } 1/ ' 2 as a function of both source 
redshift z 3 and smoothing angle 9 3 in Figs. 4-7. It follows 
the increase with redshift of the length of the line of sight. 
The variance is smaller in case of OCDM model mainly be- 
cause of the smaller normalization erg. This implies that the 
LCDM pdfs are broader compared to OCDM pdfs for a given 
smoothing angle and a given redshift. As in the case of our 
convergence studies we have a reasonable agreement from 




Figure 6. Variance plotted as a function of source redshift z 3 . 
The aperture-mass compensated filter (M ap ) is used to smooth 
the shear map. Left panel corresponds to smoothing angle 9 3 = 2' 
and the right panel corresponds to 9 3 =4'. The solid lines corre- 
spond to LCDM cosmology whereas the dashed lines correspond 
to OCDM cosmology. Data points represent averages of 10 dif- 
ferent realizations for each cosmology. Error bars are computed 
from scatter among various realizations. Dots represent LCDM 
simulations whereas OCDM simulations are represented by open 
circles. 



the fitting formula of Peacock & Dodds (1996) over the en- 
tire range of redshift and smoothing angle. Our results agree 
with previous studies which focused on the variance of com- 
pensated filters (e.g., Hoekstra et al. (2002), van Waerbeke 
et al. (2002), Jarvis et al. (2003), Hamana et al. (2003), Ben- 
abed & van Waerbeke (2003)). However, we can see that at 
large redshifts and small scales, there appears to be some dis- 
crepancy between the analytic results and the simulations. 
This deviation shows the same behaviour for the shear and 
the aperture-mass but it is larger for the latter. We must 
point out that the computation of variance only depends 
on the two-point density correlation function (or the power 
spectrum) and not on the entire correlation hierarchy or the 
parameters S p , which are the main focus of this paper. In 
other words, the deviations seen at large z 3 and small 9 3 
in Fig. 6 are solely due to the mismatch between the fit 
from Peacock & Dodds (1996) to the power spectrum and 
the numerical simulation. Therefore, they probably signal 
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9 S [in arc— minutes] fl K [in arc— minutes] 

Figure 7. Aperture-mass variance plotted as a function of 
smoothing angle 8 B . Left panel corresponds to LCDM cosmol- 
ogy whereas the right panel corresponds to OCDM cosmology. 
Three different redshifts are considered: z s = 1.5, 1 and z s = 0.5 
from top to bottom. Data points represent averages of 10 differ- 
ent realizations for each cosmology. Error bars are computed from 
scatter among various realizations. 
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Figure 8. Kurtosis (see eq. (54) for definition) of the smoothed 
shear component 7i s , as a function of the smoothing angle 9 B . 
The left panel corresponds to LCDM cosmology (solid lines) and 
the right panel to OCDM cosmology (dashed lines). Lines from 
top to bottom correspond to analytical predictions (46) (stellar 
model) for different source redshifts: z s = 0.5, 1 and 1.5. Points 
with error bars are the numerical results from simulations. 



the effects of the finite numerical resolution (note that the 
softening length increases with z). On the other hand, it 
is known that the fit from Peacock & Dodds (1996) is not 
perfect (Smith et al. 2002). 

6.3 The Kurtosis of Shear 

As pointed out in Sect. 2, for symmetry reasons all odd- 
order cumulants of shear components vanish. Therefore, the 
lowest-order non- Gaussian contribution is the kurtosis which 
we have studied as a function of both smoothing angle 6 S 
(Fig. 8) and source redshift z s (Fig. 9) (see also Takada & 
Jain 2002). We define the kurtosis as Sj is = (7? s }c/("7i s ) 3 , 
as in eq.(54). We can check that we obtain a reasonable 
agreement with the numerical simulations. There is some 
discrepancy at small scales but this might be related to the 
numerical resolution. On the other hand, it is increasingly 
difficult to predict with a good accuracy higher-order mo- 
ments of weak-lensing observables (or of the density field 
itself). In particular, while our model for the density field is 
in a sense "fitted" to the skewness S3 of the density contrast 




Figure 9. Kurtosis (see eq. (54) for definition) of the smoothed 
shear component 7i s as a function of the source redshift z s . The 
smoothing angle a is fixed at 1 arcminute in the left panel and at 
4 arcminute in the right panel. Lines correspond to the analytical 
prediction (46) (stellar model) and data points represent results 
from numerical simulations. Error bars are computed from scat- 
ter among various realizations. Although shear variance increases 
with redshift the kurtosis decreases with source redshift z a . Black 
dots with error bars represent LCDM simulation while open cir- 
cles represent OCDM simulations. 



through (25) (the quasi-linear limit S^ L is exact while the 
non-linear HEPT ansatz was seen to agree with numer- 
ical simulations), this is not the case for the kurtosis and 
higher-order moments of the density contrast. They are set 
by the simple parameterization of the generating function 
>p{y) described in Valageas, Barber & Munshi (2003) and 
Barber, Munshi & Valageas (2003) which only depends on 
5*3. Moreover, in the case of weak-lensing observables there 
is a further dependence on the angular behaviour of the cor- 
relation functions. Thus, the minimal tree-model and the 
stellar model do not give identical results for the kurtosis 
(since the former writes the four-point correlation as a sum 
of a "stellar" and a "snake" diagram while the latter only 
keeps the "stellar" graph). Nevertheless, for the shear kurto- 
sis both results are close and we only plot in Fig. 8 and Fig. 9 
the stellar model prediction since it is this model which we 
shall investigate in more details in Sect. 6.5 for the full pdf 
because it is much more convenient for numerical purposes. 

Note that whereas the variance of the shear components 
increases with the source redshift the kurtosis decreases. 
This is due both to the longer length of the line of sight 
(as we add the lensing contributions from successive mass 
sheets along the line of sight the total signal becomes closer 
to Gaussian, in agreement with the central limit theorem) 
and to the fact that the density field is closer to Gaussian at 
higher redshift. We can see that the variation with redshift 
of the kurtosis is actually very steep. This implies that any 
realistic observational study which aims to determine cos- 
mological parameters from the lower-order moments of the 
shear components will have to determine source redshifts 
very accurately. On the other hand, the analytical calcula- 
tions must take into account the spread of the distribution 
of source redshifts. This is actually straightforward within 
the formalism used in this article (see Valageas (2000a)) if 
we neglect non-linear couplings. 

In realistic surveys one must also take into account 
the finite size of the catalog as well as the observational 
noise. A complete error analysis of higher-order moments 
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Figure 10. Skewness (see eq.(54) for definition) of the aperture- 
mass Af ap plotted as a function of smoothing angle 8 S . Lines show 
the analytical prediction for z s = 0.5, 1 and z s = 1.5 from top 
to bottom (the tree-model (38) and the stellar model (46) are 
identical for the skewness). Left panel corresponds to LCDM cos- 
mology whereas the right panel corresponds to OCDM cosmology. 
Data points represent averages of 10 different realizations for each 
cosmology. Error bars are computed from scatter among various 
realizations. 




Figure 11. Skewness (see eq. (54) for definition) of the aperture- 
mass M a p plotted as a function of source redshift z s . The left 
panel corresponds to smoothing angle 9 a = 2' and the right panel 
corresponds to S =4'. The solid lines in each panel correspond to 
LCDM cosmology whereas the dashed lines correspond to OCDM 
cosmology. Data points represent averages of 10 different realiza- 
tions for each cosmology. Error bars are computed from scatter 
among various realizations. Black dots corresponds to LCDM sim- 
ulations whereas open circles correspond to OCDM simulations. 



for aperture-mass has been developed in Munshi & Coles 
(2002). Clearly from observational view points it is impor- 
tant to have a reasonable sky coverage and dense sampling of 
galaxies. Note that error analysis of lower order moments of 
shear components involves incorporating correlations among 
neighbouring cells which can be ignored for the study of 
aperture-mass statistics. This makes such studies for shear 
components more complicated although a crude order of 
magnitude study can be performed. It is also important to 
note that throughout we have ignored the noise due to in- 
trinsic ellipticity distribution which we have to include in a 
more realistic study. 

On the other hand, with future weak lensing surveys 
such as LSST and SNAP it will be feasible to study the 
non-Gaussianity induced by gravity with an unprecedented 
accuracy. It was pointed out by Hui (1999) that the skewness 
of the convergence field can directly be used to probe the 
dark energy equation of state. Future surveys with large sky 
coverage and dense sampling can also use the kurtosis of 
shear components to study the properties of dark energy. 
A clear detection of kurtosis will also help us to break the 
degeneracies in determining S7 m and as which are inherent 
in studies based on power spectrum analysis only (see also 
Takada & Jain 2002). 

6.4 Skewness of Aperture-Mass 

The shear kurtosis can be directly measured from shear 
maps hence it is a natural tool to study the departure from 
Gaussianity of the density field. However, as we noticed in 
Sect. 6.3 it suffers from several drawbacks. First, predictions 
for the kurtosis of the density field are less robust than for 
the skewness, which is a lower-order moment. Moreover, the 
relationship between the third-order moments of the density 
constrast and weak-lensing observables is stronger because 
the dependence on the angular behaviour of the correla- 
tion functions is smaller. Thus, all tree-models (including 
both the minimal tree-model and the stellar model stud- 
ied in more details in this paper) give the same predictions 



for any statistics of order three (e.g., the skewness of the 
aperture-mass or any three-point correlation) since at this 
order there is only one tree graph, whose weight is fully de- 
fined by S3. This makes the skewness of the aperture-mass 
5*3 ap the most useful probe of the non-Gaussianity of the 
density field. An additional advantage is that the aperture- 
mass provides a very localised probe of the density field in 
Fourier space, as seen in Fig. 3. This makes the compari- 
son between theory and observations more robust and more 
precise as we can probe a narrow range of scales, hence a 
well-defined regime of gravitational clustering. 

Therefore, we plot in Fig. 10 and Fig. 11 the skewness 
S 3 ap as a function of smoothing angle 9 S and source red- 
shift z s - We again obtain a reasonable agreement between 
our theoretical predictions and the numerical simulations. 
Of course, like the kurtosis Sj iB of the shear the skewness 
decreases at higher redshifts. Note however that this varia- 
tion is significantly shallower which is a useful property for 
observational purposes (the error introduced by the measure 
of the source redshifts will be smaller). 

Although present generation surveys do not provide a 
clear cosmological signal for non-zero skewness, future sur- 
veys such as Supernova Anisotropy Probe (SNAP) and the 
Large-aperture Synoptic Survey Telescope (LSST) will be 
very useful in this direction. On the other hand, various 
generalizations of the skewness S 3 ap have been proposed re- 
cently along with fast computational techniques which can 
reduce large volumes of data delivered by future weak lens- 
ing surveys. These generalisations can also handle decompo- 
sition into cosmological modes and non-gravitational modes 
induced by systematics (see e.g. Jarvis, Bernstein & Jain 
2003). 

Our studies are complimentary to other works which 
have focused on constraining or determining various halo 
model parameters from such observational studies. Our ap- 
proach is significantly different as it is based on the many- 
body correlations of the density field themselves, rather than 
relying on a decomposition of the matter distribution over 
a population of virialized halos. 
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Figure 12. Probability distribution function Vi^fis) of the 
smoothed shear component yi s . The smoothing angle 8 B is men- 
tioned in each panel while the underlying cosmological parameters 
are that of LCDM. The sources are all placed at the same redshift 
z s = 0.5. Since the pdf is even (7'(7i s ) = V{— Us)) we only display 
the pdf for positive fi s . Solid lines correspond to the analytical 
prediction (48) of the stellar model. The dashed lines show the 
Gaussian distribution with the same variance. The dotted lines 
correspond to the inclusion of the first non-zero term (kurtosis) 
in the Edgeworth expansion beyond the Gaussian approximation, 
eq.(55). Data points are results from numerical simulations. Error 
bars indicate the scatter among various realizations. 



6.5 Probability Distribution Function for Shear 
Components 7; s 

Our analytical results provide a complete analytical model- 
ing of the pdf of shear components. Extending our earlier 
studies (Valageas, Barber & Munshi 2003) we show that in- 
deed such a description is possible for the entire range of 
redshifts as well as smoothing angles of practical interest. 
Thus, we compare with numerical simulations our analyt- 
ical predictions for the pdf V(^is) of the smoothed shear 
components in Figs. 12 - 17. We consider both LCDM and 
OCDM cosmologies, for smoothing angles 8 S = 1', 2', 4' and 
8' and for three different source redshifts z s — 0.5, 1 and 1.5. 
We can see that we obtain a good agreement over the entire 
range of smoothing radius and source redshift, from linear 
to highly non-linear scales. We only plot the prediction of 
the stellar model (48) because it is much more convenient for 
practical purposes. Indeed, eq.(48) provides a simple explicit 
expression for the generating function tp liB (y)- By contrast, 
the minimal tree-model yields an implicit system (42)- (43) 
which must be solved numerically at each point y. Since the 
variable i? is actually two-dimensional this involves solving 
for a 2-d function r(i?) defined as the fixed point of eq.(43). 
This cannot be done by a simple iterative procedure for ar- 
bitary y since it would diverge for large y. Therefore, this 
implies the use of a time-consuming algorithm. Note that 
the generating function (p lis (y) must then be integrated 
over the complex plane in eq.(37). Thus, the stellar model 



Figure 13. As in previous figure but for source redshift z s = 1. 
Note that as we move towards higher redshifts the distribution 
becomes more Gaussian and the Edgeworth expansion matches 
the analytical expression better. However, truncating the expan- 
sion at any order produces spurious oscillations which are not 
physical. 
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Figure 14. As in previous figure but for source redshift z s = 1.5. 



is much simpler to implement and since both models give 
close predictions for the shear we focus on the stellar model 
in Figs. 12 - 17. Note that we have already shown in Bar- 
ber, Munshi & Valageas (2003) that both models also give 
almost identical predictions for the smoothed convergence 

We also plot the Gaussian (dashed line) and the Edge- 
worth expansion (dotted line) up to the first non-Gaussian 
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Figure 15. As in previous figure but for source redshift z s = 0.5 
and the OCDM cosmology. 



Figure 17. As in previous figure but for source redshift z s = 1.5. 




Figure 16. As in previous figure but for source redshift z s = 1. 



correction (here the kurtosis) from eq.(55). We can check 
that at large smoothing angles 9 S and source redshifts z s 
the pdf V("/i S ) becomes very close to Gaussian. This could 
already be expected from the behaviour of the kurtosis anal- 
ysed in Sect. 6.3. The departure from the Gaussian is sig- 
nificant at z s ~ 0.5 and 8 S ^ 2'. On the other hand, we 
note that the Edgeworth expansion is actually useless. It 
only provides reasonable results when the pdf is very close 
to Gaussian while as soon as there is a sizeable deviation 
from the Gaussian it introduces spurious oscillations and 
it actually fares worse than the Gaussian. Note that since 



it is an asymptotic expansion including higher-order terms 
would further worsen this discrepancy. 

One of the benefits of being able to study the shear com- 
ponents analytically is to have a means to cross check the ef- 
fects of various systematics which go into the making of con- 
vergence maps from shear data. Alternatively one can work 
with compensated filters to construct the M ap statistics as 
described above. Similarly, in a previous paper (Valageas, 
Barber & Munshi 2003) we have shown that not only the 
statistics of the shear components can be modeled in this 
manner but one can also derive the properties of the shear 
modulus |7 S |. However, we found that the correlation be- 
tween the two components is almost negligible so that the 
pdf of the shear modulus does not carry useful additional 
information. This is the reason why we have not included 
"P(|7s|) in the present study although it is a simple outcome 
of the analytical results presented here. 

Since the pdf V(^i a ) contains some information about 
the entire hierarchy of the Sp" cumulants it can be used as 
a probe of the entire pdf of the underlying mass distribu- 
tion. This is clearly apparent in eq.(48). Moreover, it may 
be used to distinguish cosmological models in a more effi- 
cient manner than relying on the few lowest-order moments. 
Given that the shear maps are a direct outcome of any weak 
lensing survey, we may hope that surveys with a low level of 
noise as well as a good sky coverage could give us some clues 
on the distribution of matter in addition to the cosmological 
parameters. However, we see in Figs. 12 - 17 that the pdf 
V(piis) is not very far from the Gaussian, partly because it 
is even so that the shape is roughly similar, hence it may be 
difficult to extract accurate measures from future surveys. 



6.6 Probability Distribution Function for 
Aperture-Mass 

As discussed in Sect. 6.4, the aperture-mass M ap shows sev- 
eral advantages over the shear components ji s . In particu- 
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Figure 18. Probability distribution function 'P(M a p) for the 
Map statistic which uses a compensated filter. We probe vari- 
ous smoothing radii from 9 S = 2' up to 9 a = 16'. All sources are 
placed at source redshift z B = 0.5 and we consider the LCDM cos- 
mology. Solid lines show the prediction of the stellar model (48). 
Long dashed lines correspond to the minimal tree-model (42)- 
(43). Short dashed lines are the Gaussian with the same variance 
while the dotted lines are the Edgeworth expansion (55) includ- 
ing the first correction to Gaussianity (skewness). Data points 
are results from numerical simulations. Ten different realizations 
are analysed and error bars denote the scatter among various re- 
alizations. Note that for negative M ap the minimal tree-model 
predictions are much closer to numerical simulations than the 
stellar model. 
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Figure 20. As in previous figure but for source redshift z s = 1.5. 
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Figure 19. As in previous figure but for source redshift z s = 1. 
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Figure 21. As in previous figure but for source redshift z s =0.5 
and the OCDM cosmology. 



lar, it probes a narrower range of scales and it is not an even 
quantity. Therefore, we plot in Figs. 18 - 23 the pdf 'P(Afap) 
of the aperture-mass. We again consider both LCDM and 
OCDM cosmologies and sources redshifts z s = 0.5, 1 and 1.5. 
However, we now study the smoothing angles 8 S = 2', 4', 8' 
and 16'. Indeed, as shown in Fig. 3 the aperture-mass probes 
smaller scales than the shear for a same angular radius, 
hence the limitation from the finite resolution of numeri- 
cal simulations shifts to larger angles. On the other hand, 
we now plot the analytical predictions of both the mini- 
mal tree-model (42)- (43) (long dashed lines) and the stellar 
model (48) (solid lines). Indeed, contrary to the convergence 
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Figure 22. As in previous figure but for source redshift z s = 1. 
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Figure 23. As in previous figure but for source redshift z s = 1.5. 



(see Barber, Munshi & Valageas 2003) or the shear, the pre- 
dictions of both models significantly differ at smaller angu- 
lar scales for the large negative M ap tail. As explained in 
Sect. 4.4, this behaviour can be understood on simple the- 
oretical grounds. At very small angles (9 S ^ 0.4') the pdf 
P(M ap ) predicted by the stellar model shows a singular be- 
haviour and it vanishes for Af ap < 0. At the scales of pratical 
interest displayed in Figs. 18 - 23 the pdf is still regular but 
we can see a hint of the trend towards this singularity as 
the falloff at negative M ap becomes increasingly steep for 
smaller scales. By contrast, the pdf obtained from the mini- 
mal tree-model which does not suffer from this singular be- 



haviour shows a smoother cutoff at negative M ap (although 
it remains sharper than for positive M ap , in agreement with 
numerical simulations) . Note that this behaviour shows that 
no realizable density field can exactly obey the stellar model 
(whatever the coefficients S p ), while this remains an open 
question for the minimal tree-model. 

We can see in the figures that the prediction from the 
minimal tree-model shows a good agreement with the nu- 
merical simulations over all scales and redshifts of interest. 
On the other hand, as expected from the previous discus- 
sion the stellar model fails to reproduce the tail at negative 
M ap , except for large angles and redshifts. However, it yields 
a reasonable prediction for the pdf for M ap > (albeit the 
minimal tree-model fares slightly better in this domain too). 
Therefore, we can conclude that the minimal tree-model is 
a much better tool to study the statistics of the aperture- 
mass. Fortunately, although the numerical computation is 
still more difficult and more computer time consuming than 
for the stellar model, the function r(i?) = t(i?) introduced in 
eq.(43) is now only one-dimensional which makes the com- 
putation much easier than for the shear (where t(i?) was 
truly 2-dimensional) , since the filter J/M ap ($) obeys a radial 
symmetry. 

We also display the Gaussian (dashed lines) with the 
same variance. Of course, at large redshifts and angles the 
pdf becomes closer to Gaussian but we can see that even for 
9 a ~ 8' and z s ~ 1.5 the deviations from Gaussianity are 
clear. Moreover, at smaller angles and redshifts the Gaus- 
sian completely fails to reproduce the pdf obtained from 
numerical simulations. It cannot follow the sharp peak at 
M ap ~ 0, the extended tails and the asymmetry of the pdf. 
Thus, we can see that the departure from Gaussianity is 
much more important for "P(M ap ) than for the pdf V{"fia) 
of the shear studied in Sect. 6.5. Therefore, the pdf P(Map) 
of the aperture-mass should provide a much more efficient 
tool to measure the deviations from Gaussianity than the 
pdf V(~iis) of the shear. Note that non-Gaussian features, 
like low-order moments for instance, exhibit a strong fl m - 
dependence, mainly due to the presence of the normalising 
factor in eq.(l) (e.g., Bernardeau et al. 1997). Therefore, the 
aperture-mass is a convenient tool to measure such cosmo- 
logical parameters as well as the properties of the underlying 
density field. On the other hand, we also plot the Edgeworth 
expansion (dotted lines) up to the first non-Gaussian term 
(skewness). As was the case for the shear, we can see that 
this asymptotic expansion is of very little use, since it only 
fares well at large angles and redshifts where the Gaussian 
is already a reasonable description (although it somewhat 
improves the shape of the pdf near its maximum) while as 
soon as there is a significant deviation from the Gaussian it 
yields spurious oscillations which make it useless (in some 
domains it even gives negative values for P(M ap )). 

Attempts in modeling the pdf V(M a , p ) were initiated 
by Reblinsky et al. (1999), where a halo model of clustering 
was assumed in order to compute the positive tail of the 
pdf. Indeed, the high-M ap tail is very efficient in probing 
and mapping out large concentrations of dark mass in weak 
lensing surveys. On the other hand, a complete prediction of 
the pdf for all values of M ap , which would be difficult to ob- 
tain from a halo model which cannot faithfully describe low- 
density regions and filamentary structures, gives a valuable 
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Figure 24. Probability distribution functions "Pfais) of the 
smoothed shear component -yi s (right panel) and r P{M ap ) for the 
Af a p (left panel) are compared for two different cosmologies for a 
given smoothing angle 9 S = 2' and source redshift z B = 1. Solid 
lines correspond to LCDM cosmology and dashed lines correspond 
to OCDM cosmology. 

insight in probing both overdense and underdense regions 
in a statistical manner. We have shown that our method, 
which is based on the many-body correlations rather than 
on a decomposition over halos, provides such a model. In 
particular, the results presented above show that the neg- 
ative M a p-tail is very sensitive to the angular behaviour of 
the correlation functions (in addition to their amplitude). 
Therefore, it provides a direct probe of higher order cor- 
relations and it can help us to discriminate among various 
models of gravitational clustering. In fact, we have shown 
that it already rules out the stellar model even though this 
model provides a very good description for the statistics of 
the convergence (Barber, Munshi & Valageas 2003) and the 
smoothed shear (Sect. 6.5). 

Our results extend the work by Bernardeau & Valageas 
(2000) who presented a detailed study of the pdf P(M ap ) 
within both the quasi-linear regime (where exact calcula- 
tions are possible) and the highly non-linear regime (where 
they used the minimal tree-model). They also compared the 
predictions of the minimal tree-model with numerical sim- 
ulations for various cosmologies, with z s = 1 and 9 a =4'. 
Here we have shown that the minimal tree-model actually 
provides good predictions for P(M ap ) over all scales and 
redshifts of interest, from quasi-linear to highly non- linear 
scales. Moreover, by introducing the stellar model we have 
shown that V(M ap ) could be used to probe the detailed an- 
gular behaviour of the many-body correlations in addition 
to their amplitude. A more complete study of "P(M ap ) with 
realistic source distributions as well as noise due to the in- 
trinsic ellipticity distribution will be presented elsewhere. 

6.7 Comparison of the Probability Distribution 
Functions obtained for both cosmologies 

Finally, for the sake of completeness, we plot in Fig. 24 
the pdfs obtained in both cosmological models for the shear 
and the aperture-mass. The difference between both models 
is rather small. Indeed, these weak-lensing observables are 
mostly sensitive to the matter density parameter Q m and 
to the amplitude as of the density fluctuations while the 
dependence on Qa is smaller. However, note that the com- 
parison is not so straightforward since both Qa and as differ 



between the two models (and most of the change of the pdf 
comes from as)- A more precise study of the dependence on 
cosmological parameters and on the source redshift is left 
for future works. 



7 DISCUSSION 

Weak lensing surveys are being regularly used to probe the 
matter distribution and the underlying cosmology. However, 
the modeling of statistical quantities related to weak lensing 
surveys, such as convergence maps or related shear maps, is 
rather difficult because of the non-linear nature of the un- 
derlying density field they try to probe. Besides, the non- 
linearity grows at smaller angular scales which should dom- 
inate the signal in forthcoming surveys. In a series of re- 
cent papers we have shown how to combine a widely used 
parameterization of the non-linear matter power-spectrum 
with various simple models of the matter density correlation 
hierarchy. This allows a detailed description of weak lensing 
statistics. 

In the present paper we mainly focus on the statistical 
properties of the weak lensing shear components 7i B and the 
closely related aperture-mass M ap . These two quantities can 
be computed from observational surveys in a direct manner 
from galaxy ellipticities (in the weak lensing regime). By 
contrast, extracting convergence maps is more difficult as it 
involves a non-local inversion problem and exhibits a mass 
sheet degeneracy. However, we must note that the measure 
of shear statistics still remains a difficult task because of the 
non-trivial survey topology (a fraction of the survey area has 
to be excluded from the analysis because of image defects, 
bright stars,..). Thus, rather than computing the shear or the 
aperture-mass realized over many complete circles on the 
sky, one usually estimates low-order correlation functions 
of the shear which are next integrated in order to obtain 
low-order moments of the smoothed shear or aperture-mass. 
However, we do not study these points in this article. We 
investigate two simple analytical descriptions of the matter 
density field: the minimal tree model and the stellar model 
which we introduced in Valageas, Barber & Munshi (2004). 
Both models give identical results for the statistics of the 
3-d density contrast smoothed over spherical cells and only 
differ by the detailed angular dependence of the many-body 
density correlations. 

In a previous study (Barber, Munshi & Valageas 2004) 
we have shown that both models also give almost identi- 
cal results for the smoothed convergence k s . In agreement 
with earlier works, this shows that the pdf V(k s ) (or its mo- 
ments) provides a robust measure of the pdf V(Sr). This 
also enables one to measure in a robust fashion the un- 
derlying cosmological parameters as well as the amplitude 
of the density correlations. In the present study we extend 
such calculations to more intricate compensated filters which 
can directly be constructed from shear maps: the smoothed 
shear components 7^ and the aperture- mass M ap . As these 
observables involve more intricate filters we can expect the 
dependence on the angular behaviour of the many-body cor- 
relations to be larger than for the smoothed convergence. We 
find that both models actually yield rather close results for 
the smoothed shear components and the positive tail of the 
aperture-mass. Moreover, we obtain a good agreement with 
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numerical simulations over all scales and redshifts of prac- 
tical interest. Therefore, in this domain the shear and the 
aperture-mass provide again a robust constraint on cosmo- 
logical parameters and the statistics of the smoothed 3-d 
density contrast Sr. Besides, we found that V(M ap ) shows a 
stronger departure from the Gaussian than P(jis), so that 
the aperture-mass appears to be a very useful tool in this 
respect. 

On the other hand, we also note that at small angles 
(9 S 5s 2') the tail of the pdf V(M ap ) for negative M ap shows 
a strong variation between both models. Whereas the mini- 
mal tree-model provides a good description of the numerical 
data down to the smallest scales available to us the stellar 
model shows a significant discrepancy (while the part of the 
pdf over M ap > remains reasonable). The stellar model 
actually breaks down at 9 S <> 0.4' for M ap < 0. This clearly 
indicates that the aperture-mass statistics can be used as 
a very precise probe of the high-order correlation hierarchy. 
Contrary to the smoothed convergence k 3 , it is not only sen- 
sitive to the amplitude of the many-body correlations but 
also to their detailed angular behaviour. Thus, it provides 
a complimentary tool to the smoothed convergence or shear 
components. From a theoretical point of view, we can note 
that this behaviour also means that no physical density field 
can be exactly described by the stellar model (while this 
remains an open issue for generic minimal tree-models) al- 
though it provides a very interesting and convenient approx- 
imation which works very well for both the smoothed con- 
vergence and the smoothed shear components. This shows 
that the shape of the pdf "P(M ap ) can actually rule out sen- 
sible models of gravitational clustering. 

Combined with our previous results for shear and con- 
vergence statistics (Valageas, Barber & Munshi 2004, Bar- 
ber, Munshi & Valageas 2004) we have developed a very 
powerful technique to analyse weak lensing survey results. 
Although in our present calculations (both analytical and 
numerical) we have not included observational details like 
the finite width in source distribution or the noise due to the 
intrinsic ellipticity distribution of galaxies, they can be in- 
corporated easily in our computations. A detailed and more 
elaborate comparison will be presented elsewhere when such 
simulations become available. 

Several interesting numerical issues became clear from 
our study. For a given angular scale 9 S compensated filters 
like UM av pick up more contributions from higher wavenum- 
bers as compared with tophat filters. This makes them more 
difficult to study numerically at smaller angular scales, as 
the finite resolution of the simulations starts to play a role. 
This is also the regime where M ap statistics are very sen- 
sitive to the detailed analytical modeling of the correlation 
hierarchy. Our work serves as a precurser to studies using 
simulations with much finer resolutions. However, using such 
compensated filters also means that finite volume effects are 
much less pronounced as compared with tophat filters (see 
Munshi & Coles (2002) for a detailed analysis of various 
spurious results in determination of M ap statistics). 

Most analytical studies in the non-linear regime use a 
halo model (for detailed predictions of a halo model see 
Takada & Jain 2003). In principles, halo models (see Cooray 
& Sheth 200 for a review) can predict the higher-order corre- 
lation functions in the highly non-linear regime. Indeed, on 
small scales the latter are set by the density profile of the 



halos as the p— point correlations are dominated by the con- 
tribution associated with all p points being within the same 
halo. However, it is interesting to note that the neglect of 
substructures may lead to larger inaccuracies for high-order 
statistics. On the other hand, at intermediate scales one also 
probes the correlations among different halos which intro- 
duces new unknowns and makes explicit calculations cum- 
bersome. This is important for handling projection effects 
which involve the mixing of various scales. Finally, such a 
model for the density field is not well-suited to describe the 
low-density and underdense regions (e.g., voids, filaments) 
which are outside virialized objects. Our method follows a 
completely different approach based on the high-order corre- 
lation functions themselves rather than on a decomposition 
of the matter distribution over a population of virialized ha- 
los. Clearly such an approach provides an independent com- 
plimentary scheme and will be useful for various numerical 
cross-checks. 

The simulation technique that we employ is completely 
different from the more popular ray tracing methods. Our 
studies therefore not only provide a good test for analytical 
results but it is also a good consistency test for such new 
simulation techniques. 
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